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Abstract 

We examine the structure of infrared singularities in QCD and quantum General Relativity, 
from the point of view of the recently conjectured double copy property which relates scattering 
amplitudes in non-Abelian gauge theories with gravitational counterparts. We show that IR 
divergences in both theories are consistent with the double copy procedure, to all orders in per- 
turbation theory, thus providing all loop-level evidence for the conjecture. We further comment 
on the relevance, or otherwise, to the so-called dipole formula, a conjecture for the complete 
structure of IR singularities in QCD. 



1 Introduction 

Scattering amplitudes in quantum field theories continue to generate a large amount of interest, 
due to the many phenomenological and theoretical applications. Much of this work (motivated by 
the relevance of QCD to hadron colliders) focuses on non-Abelian gauge theories, including super- 
symmetric counterparts (e.g. A/" = 4 Super- Yang Mills theory) whose simpler structure provides a 
useful testing ground for theoretical techniques. Amplitudes have also been studied in gravitational 
field theories, namely General Relativity and its supersymmetric counterparts such as = 8 super- 
gravity, which may provide an ultraviolet finite field theory of gravity [Tj. Although gravitational 
theories look superficially very different from non-Abelian gauge theories, a recent programme of 
work has suggested that gauge and gravity amplitudes may be related in an intriguing way [2j|4]. 
Firstly, gauge theory amplitudes may be written in a special form which manifests an explicit 
duality between colour and kinematics, the so-called BCJ duality of |5]. Secondly, an m-point, 
L-loop order gauge theory amplitude in BCJ dual form can be translated into an equivalent gravity 
amplitude by the double copy procedure, in which colour factors are replaced by kinematic factors. 
That gauge theory amplitudes admit BCJ duality is known to be the case at tree level (3|[6[|9], 
where the double copy property is equivalent to the known KLT relations (To| relating amplitudes 
in gauge and gravity theories. At loop level there is no formal proof of the double copy conjecture, 
although it is known that this should hold pending the existence of BCJ duality to all orders in the 
gauge theory, for the cases of pure gravity and Af = 8 super-gravity [s] . Loop-level checks of the 
duality have been carried out up to four loop order in A/" = 4 SYM theory [2 , 11 - 13 for amplitudes 
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involving up to five external particles, and two loop order (for four point scattering with maximal 
helicity violation) in QCD |2]. 



As is perhaps clear from the above comments, much work on BCJ duality and the double copy 
property has focussed on = 4 SYM theory, in which a restricted set of loop diagrams con- 
tribute f4l . A significant motivation for studying AA = 4 SYM is that the double copy relates this to 
M = 8 supergravity, so that one may investigate the issue of whether the latter theory is ultraviolet 
finite [T 14 ,15 . To date, less work has been invested in examining the non-supersymmetric context 
of pure Yang-Mills theory, where the double copy relates this to General Relativity. Examples 
include 16,17] (see also |18|), in which amplitudes in 4 < AA < 8 supergravity are considered by 
combining BCJ dual AA = 4 SYM amplitudes with those in a less supersymmetric theory, exploiting 
the fact that the two sets of kinematic factors in the double copy procedure need not both satisfy 
BCJ duality. Recently, an extension of the duality was considered, to Yang-Mills theory deformed 
by higher dimensional operators 
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The aim of this paper is to point out that the infrared limit of QCD [^provides a clean environment 
in which to examine the double copy property in a non-supersymmetric context. We will argue that 
the known structure of infrared singularities in both QCD and gravity is consistent with the double 
copy, via BCJ duality, to all loop orders. The infrared behaviour of gauge theory amplitudes has 
been intensively studied over many years (see e.g. [20f[29] ) . In QCD, the state of current knowledge 



regarding massless scattering amplitudes can be summarised in the so-called dipole formula 130 - 34 



a conjecture which states that soft and collinear singularities exponentiate such that the exponent 
contains colour correlations between at most pairs of external legs. The formula is known to be 
exact up to two loops, and the structure of possible corrections at three loop order and beyond has 
been further studied in [35 -39 . It is known the dipole formula fails already at two loop level for 
massive external particles 40-48 , a fact which is not immediately relevant here due to the fact that 



the framework of BCJ duality and the double copy procedure is set up only for massless particles. 



The study of infrared singularities in General Relativity was first examined in 49 , and there has 



recently been a revival of interest, particularly in the work of 50-52 which aims to describe soft 
graviton physics in the same language used in a gauge theory context, and which we will review 
in what follows. It is now well established that the soft limit of gravity is one loop exact. That 
is, infrared singularities exponentiate such that the exponent contains only one loop graphs. Note 
that this implies that all IR singularities in gravity are dipoles, in the sense that at one loop in the 
exponent only pairs of particles can be correlated. 



The above comments, and the relationship between gauge theory and gravity offered by the double 
copy procedure, beg the question: could the QCD dipole formula have an essentially gravitational 
origin? This thought motivated the present study, and we will see that in fact the answer appears 
to be no. Nevertheless, we will find that the structure of IR singularities in QCD match up at 
all orders with the known structure of singularities in General Relativity, providing all-loop-level 
evidence for the double copy conjecture. The argument is insensitive to possible corrections to 
the dipole formula - we will see explicitly that many singularities disappear when the double copy 

^We consider pure gluodynamics (QCD without quarks) throughout the paper. For convenience, we refer to this 
as QCD throughout. 
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is taken. This includes collinear singularities, which are already known to vanish in gravity after 
summing over diagrams 49 , 52 . 



The structure of the paper is as follows. In section [2] we review necessary background material in 
more detail, namely concepts relating to BCJ duality, the double copy, and the structure of infrared 
singularities in QCD and gravity. In section [3] we examine BCJ duality and the double copy in 
the soft limit at one loop order, before extending this analysis to two loop order in section [4j In 
section[5]we generalise our remarks to all loop orders, before summarising our results and concluding 
in section [6l 



2 Review of necessary concepts 
2.1 BCJ duality and the double copy 

In this section we briefly summarise BCJ duality [5j and the associated double copy conjecture (2][3], 
which together provide a map from gauge theory to gravity amplitudes, potentially at the multiloop 
level B 

A general massless m-point L-loop gauge theory amplitude Am^ in D space-time dimensions can 
be written as 

a(L) _ -L m-2+2LSr^ [ Y] ^^Pl 1 rijCj 



where g is the coupling constant. Here the sum is over the complete set of graphs involving triple 
gluon vertices, consistent with the given loop order and number of external particles, and Si a 
symmetry factor for each graph i (the dimension of its automorphism group). The denominator 
contains all relevant propagator momenta, and rij is the kinematic numerator associated with each 
graph. Finally, Cj is the colour factor of each graph, obtained by dressing each three gluon vertex 
with a factor (adopting the same conventions as [2]) 

where f"'^'^ are the usual SU(3) structure constants. Noteworthy points regarding this formula are: 

• The restriction to graphs with cubic vertices only does not mean that graphs with four gluon 
vertices have been excluded. Rather, one can always choose to rewrite the latter in terms 
of cubic graphs, by introducing extra momenta in the denominator (which are compensated 
by additional factors in the numerator). This relies upon the fact that the colour factor for 
a four gluon vertex is a product of two structure constants, consistent with a pair of three 
gluon vertices. 

• The numerators in eq. ([T]) may or may not come from individual Feynman diagrams. They 
may also have been obtained from e.g. a generalised unitary-based approach (see e.g. [53] 
and references therein), in which case each rij collects all kinematic information associated 
with a particular scalar integral (fixed by the denominator). 

^See also section 11(A) of [4j for a recent and pedagogical review of this material. 



3 



Figure 1: Illustration of three graphs related by BCJ transformations, where the grey boxes denote 
the rest of the amplitude. 



From any graph at a given loop order, one may construct two other graphs related by taking BCJ 
transformations involving crossing of s-channel like subgraphs into t and u-channel subgraphs. 
This is best illustrated pictorially, as in figure [T] Such transformations exist for any internal line 
(propagator) in a given cubic graph, leading to a set of inter-related graphs, which can be generated 
from a single diagram. The colour factors associated with the three graphs in figure [T] satisfy the 
identity 

Cs = Ct + Cu, (3) 

where Cj is the colour factor associated with the diagram containing an i-channel-like subgraph. 
This follows from the Jacobi identity for the structure constants, after factoring out the part of 
each colour factor which is the same for each graph. Note there is an ambiguity in how one defines 
signs in this identity (which can be compensated by a corresponding choice for the numerators). 
We here use the same conventions as [4]. 

It is conjectured that one can always choose to redefine the numerators nj in the amplitude of 
eq. ([l]), such that they satisfy a similar identity to eq. ([S]). That is, the numerators of the s, t and 
ti-channel-like graphs in figure [T] can be chosen to obey 

ns=nt + nu. (4) 

Furthermore, if a given colour factor picks up a minus sign under interchange of two (external or 
internal) legs, then 

Ci -Ci nj -)> -rij. (5) 

This properties are collectively known as BCJ duality, after |5|. It has been proven at tree-level, but 
remains a conjecture at loop level. Transformations required to write the numerators in a suitable 
form are known as generalised gauge transformations, and an algorithmic procedure exists in prin- 
ciple to establish, for a given loop-level amplitude, what a possible set of numerators actually is [ij. 
It is not clear, however, whether this algorithm is fully general, particularly in non-supersymmetric 
theories where more diagrams enter than in supersymmetric cases. In general, one may represent 
the effect of a generalised gauge transformation on a given numerator Ui as 

nj-^n-i + Aj, (6) 

where the {Aj} must satisfy 



E 



=0. (7) 
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Again the sum is over all cubic diagrams with colour factors Cj and propagator denominators . . 
Indeed, eq. Q is obtained by substituting the transformation of eq. ^ into eq. 

and requirmg 

that the amplitude is invariant. One may further decompose the generalised gauge parameters Aj 
(assuming these to be local) as [3] 

Ai = Y,^i,a,Pl,- (8) 

That is, the Aj associated with diagram i can be expanded in terms of the inverse propagators of 
this diagram. The Aj factors then move contributions between diagrams by cancelling propagators, 
so that contributions from one diagram can be absorbed in another. 



Once a form of the amplitude satisfying BCJ duality has been found, a remarkable conjecture states 
that one can turn it into a gravity amplitude, in a straightforward manner. By stripping off the 
colour factors in eq. and replacing them with another set of numerators {hi}, the double copy 
conjecture states that |2 , 3p" 



is an m-point, L-loop gravity amplitude. Note that the numerators {ui} and {hi} need not come 
from the same theory. In this paper, we will be concerned with both sets of numerators coming 
from QCD, in which case the gravity amplitude corresponds to general relativity (coupled to an 
antisymmetric tensor and a dilaton). If the two gauge theories are M = N and N = M Super- 
Yang-Mills theory, then one obtains an amplitude in A/" = {N + M) supergravity. 

This is all we need for what follows. We now turn to the study of infrared singularities in QCD 
and gravity. 



2.2 Infrared singularities in non-Abelian gauge theory 

Infrared singularities in quantum field theory have been studied over many decades, with a vast 
accompanying literature. Here we briefiy summarise only those facts which are of direct relevance 
to what follows. For a more detailed pedagogical review, see e.g. |34j. Whilst our statements can 
be interpreted in the context of a general non-Abelian gauge theory, we explicitly refer to QCD 
throughout. 

It is by now well-known that infrared singularities factorise, such that the general structure of an 
m-point scattering amplitude in QCD is, schematically, 

m J 

Am='Hm-S-J{^. (10) 

, Ui 
1=1 

Here Hm is an infrared finite hard interaction, which is dressed by a universal soft function S that 
collects all infrared singularities. The jet function Ji collects collinear singularities associated with 



*Note that in this paper we use k = ^WttGn rather than k = \/327rG]v (see section 2.31. This modifies the 
couphng factors in eq. ([9| relative to those in pllsl. 
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external leg i, and the eikonal jet function Ji removes the double counting of divergences which are 
both soft and collinear. The soft function is given by a vacuum expectation value of Wilson line 
operators along the space-time trajectories of the outgoing hard particles. Equivalently, one may 
calculate the soft part of a given hard interaction by dressing all external lines with all possible 
soft gluon emissions, according to the eikonal Feynman rule 



9s T, 



p ■ 



(11) 



where gs is the strong coupling constant, and p (k) the momentum of the hard external line (soft 
gluon) respectively. Furthermore, Tj is a colour generator in the representation of external line i, 
where we have used the notation of 54,55 . Where soft gluons meet off the external lines, they 



couple according to the usual three and four gluon vertices of QCD. 



All singularities appearing in eq. ( 10 ) can be shown to exponentiate, so that instead of eq. ( 10 ) we 
may write 

Am = T~Lm • Z, (12) 



where 



exp 



^Cn{{pi],e,n)a'^s 



n=l 



(13) 



contains all soft and collinear singularities, and depends upon the momenta {pj}, as well as the 
dimensional regularisation parameter e and scale /i. The structure of the exponent (i.e. the form 



of the coefficients {cn} in eq. (13)) is known explicitly up to two loop order for both massless and 



massive particles [40 -48 ,[56, 57 . For massless particles, it has the remarkable property of involving 
both colour and kinematic correlations between at most pairs of particles, despite the fact that one 
would naively expect correlations between more than two particles to appear at two loop order and 
beyond. This property motivated the conjecture of the so-called dipole formula in QCD pO-34 



which gives the all-order structure of eq. (13) as 



exp < 



dA2 
A2 



A2 



^ m 



7 J, (a5(A^e)) 



(14) 

Here jk is the cusp anomalous dimension (itself a perturbative expansion in as with constant 
coefficients), and 7j. a further anomalous dimension associated with jet i, and which collects hard 
collinear contributions. The double sum in the first term is over all pairs of particles {i,j), and 
following 



34 



we have used the notation Ajj = 1 ii i and j are both in the initial or both in the final 



state (otherwise Ay = 0). Equation (14) indeed contains correlations between dipoles only, and is 
known to break down already at two loop order for massive external legs. For the massless case, 
corrections may occur at three loop order and beyond, either through explicit kinematic dependence 
of the relevant Feynman integrals, or through higher order Casimir invariants appearing in the cusp 
anomalous dimension. The form of possible corrections has been investigated in 35 -39 . Further 



progress may be made using recently developed techniques for classifying the structure of the 



exponent 28,29,58,59 , or by considering alternative gauges [60j . 



6 



2.3 Infrared singularities in gravity 

Complementary to the above mentioned studies in gauge theory, IR singularities have also been 
investigated in gravity, commencing with the classic work of |49| . Recently, there has been a revival 
of interest, which has focused in particular on writing the structure of gravitational IR divergences 
using the same language as is used in modern QCD studies. Reference \50i suggested the use of 



the following gravitational generalisation of eq. (10): 



Mm = -Hm-S^'''''-, (15) 

where A4m is an m-point gravity amplitude. Here Tim is again a hard interaction which is infrared 
finite, and S^'^'^^' is a universal gravitational soft function, which is given by a vacuum expectation 
value of suitable Wilson line operators. There are no jet functions, due to the fact that collinear 
singularities cancel in gravity after summing over all diagrams and using momentum conservation. 
The latter property was first established in the soft limit [49] , but has recently been fully ex- 
tended to encompass hard collinear singularities 52,61 . Furthermore, ref. [511 examined the form 



of eq. (15) in more detail using the path integral approach of |62j , also classifying what happens 



beyond the eikonal approximation. A similar structure of next-to-eikonal corrections was found as 



in the case of gauge theory, as explored in 62 , 63 . Gravitational Wilson lines in a soft-graviton 



context were further explored in [64] using the radial coordinate space picture of [60l , and a continu- 
ous deformation was found between the cusp anomalous dimensions of QED and gravity at one loop. 

As in the gauge theory case, the soft function in gravity exponentiates. However, a drastic simplica- 
tion over gauge theory occurs in that the exponent contains only one-loop diagrams, with no higher 
order corrections. This property is known as one-loop exactness, and has been firmly established by 



the studies of (49 50 ,52 . It implies that all infrared singularities (i.e. to all orders in perturbation 



theory) ultimately stem from the exponentiation of the one-loop corrections, in marked contrast to 



the gauge theory case of eq. (13): even if the dipole formula of eq. (14) happens to be true, there 
is still a further perturbation expansion to be carried out in the exponent, whose soft part requires 
the cusp anomalous dimension to all orders. It is amusing to note, as has already been mentioned 
in the introduction, that one-loop exactness in gravity implies that all infrared singularities are 
associated with pairs of particles (the most that can be correlated with a single graviton exchange) . 
This is reminiscent of the QCD dipole formula in some sense (i.e. that higher multipole correlations 
vanish), and allows us to speculate as to whether the dipole formula may have a gravitational origin. 
The results of this paper would appear to suggest that this is not the case, due to the disappearance 
of many singularities upon taking the double copy of a gauge theory. However, one-loop exactness 
has another important role: it tells us that we know the all-order structure of IR singularities in 
gravity completely. We can then ask whether the known IR divergence structures in QCD and 
gravity are consistent with each other, if we apply the double copy procedure (via BCJ duality). 
If this is so, this provides all-order evidence (at least in a particular limit) for the double copy 
conjecture. This will be the aim of the rest of the paper. 

In our subsequent calculations we will use the following conventions (a number of different choices 
exist in the literature - see e.g. [65| for a convenient reference). The graviton field h^i, is defined in 
terms of the metric tensor g^^i, via 

9iiu = rjiiu + Khfj^y, (16) 
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where rj^i, is the Minkowski space metric, k = \/167tGn and is Newton's constant. Emission 
of soft gravitons with momentum k from a hard external hne of momentum p is then given by the 
eikonal Feynman rule (see 64 for a recent derivation using the above conventions) 



2 p-k 



(17) 



which may be compared with the gauge theory eikonal Feynman rule of eq. (11). Finally, we will 
need the graviton propagator, for which we use the result in the de Bonder gauge: 



-iP, 



A:2 



D 



(18) 



Having reviewed the necessary material for what follows, we begin our investigation of the double 
copy procedure in the soft limit in the following section. 



3 One loop analysis 

In the previous section, we reviewed various theoretical ideas concerning BCJ duality, the double 
copy, and the structure of infrared singularities in both QCD and gravity. This motivated the 
central question of our study: do the known IR singularity structures in QCD and gravity match 
up with each other upon taking the double copy of QCD? The aim of this paper is to argue that 
this is indeed the case, and we will begin at one-loop order with the following strategy: 

1. We classify all possible BCJ relations of the form of eq. Q, between sets of three diagrams 
containing cubic vertices, in the case of pure Yang-Mills theory. 

2. Next, we write down all possible diagrams which give infrared singularities at this order, and 
show that they can be matched up in sets of three, each of which offers an explicit solution 
of a BCJ relation in the soft limit. 

3. Having shown that BCJ duality is satisfied (up to corrections which are subleading in the 
soft limit), we then take the double copy of all infrared divergent diagrams, and verify that 
this reproduces the one-loop divergences in GR. 

The first step is to obtain the BCJ relations, which we do in the following subsection. 
3.1 BCJ relations 

To start with, we must write down all possible one-loop graphs, and for a concrete example we 
consider the case of 4-point scattering, labelled as in figure [2| We may also define the usual 
Mandelstam invariants 

s = {pi+P2f, t = {pi-Vzf, u = {pi-pif. (19) 
The box and triangle graphs are labelled in figure [3j The bubble graphs are shown in figure |4j 

As explained in section [2| the diagrams in figures [3] and |4] form sets of three, related by subjecting 
internal lines to a BCJ transformation (replacement of t-channel-like exchange by s- and n— channel 
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Figure 2: Four point scattering, where all momenta are taken to be outgoing. 




Figure 3: Set of all box and triangle graphs at one-loop order, where the relevance of the numeric 
labels is explained in the text. 
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(b1) (b2) (b3) (b4) 




(b13) (b14) (b15) 



Figure 4: Set of bubble graphs at one-loop order. 



exchanges). The colour factors of such a set satisfy the Jacobi identity, and BCJ duality, if satisfied, 
then requires a corresponding equation to hold for the kinematic numerators associated with each 
graph. The BCJ relations in the present case can be split into two kinds. Firstly, there are relations 
which relate two box topologies to a triangle. These are (using the labels from figures [3] and [4]) 



-UBi + nB2 + riTi 


= 


-UBl + nB2 + nT3 


= 


-UBI + nB2 + nT5 


= 


-UBl + nB2 + nT7 


= 


-TLBS + nB2 + nT9 


= 


-nss + nB2 + riTii 


= 



(20) 



Next, there are BCJ relations that relate two triangle topologies with a bubble. There are two 
subclasses - firstly those relations which involve an external self-energy: 





;i 


2,3,4,5) + nTi(l,3,2,4,5) 


- ribi 


= 


riTi 


;i 


2,3,4,5) + nTi(l,2,5,4,3) 


- nb2 


= 


nT2 


;i 


2,3,4,5) + nT2(l,3,2,4,5) 


- 


= 


nT2 


;i 


2,3,4,5) + nT2(l,2,5,4,3) 


- nbi 


= 


UTS 


;i 


2,3,4,5) + nT3(l,3,2,4,5) 


- ribe 


= 


riTs 


;i 


2,3,4,5) + nT3(l,2,5,4,3) 


- ribs 


= 



10 



nTzifl 2 3 4 5) + nTdd 3 2 4 5)- 




= 


'nm.(^ 9 4 i^"; -I- Tir^- n 9 1^ 4 '^'^ 
'1T4VJ^) ^) <5) 4) Oj + ''T4(,J-) ^) 0) ^; "Jj 




— n 


nT5(l,2,3,4,5) + nT5(l,3,2,4,5)- 


mil 


= 


nr5(l,2,3,4,5) + nT5(l,2,5,4,3)- 


nbi3 


= 


nT6(l,2,3,4,5) + nT6(l,3,2,4,5)- 




= 


nr6(l,2,3,4,5) + nT6(l,2,5,4,3)- 




= 



(21) 

Here we have been careful to label the momenta associated with each leg in the triangles, as shown 
in figure [3} Effecting a BCJ transformation on internal lines in the triangle graphs can "twist" 
the triangle. The resulting numerator factor may be written in terms of the numerator of the 
original triangle topology, but with leg labels interchanged. For brevity, we do not bother labelling 
momentum legs for the bubble and box graphs. 

Secondly, there are those which involve an internal self-energy: 



nTi(l,2,3,4,5) + nTi(l,5,3,4,2] 


- nfes 


= 


nT2(l,2,3,4,5) + nT2(l,5,3,4,2) 


- nfes 


= 


nT3(l,2,3,4,5) + nT3(l,5,3,4,2) 


- "-bio 


= 


nT4(l,2,3,4,5) + nT4(l,5,3,4,2) 


- "-wo 


= 


nT5(l,2,3,4,5) + nT5(l,5,3,4,2) 


- nbi5 


= 


nT6(l,2,3,4,5) + nT6(l,5,3,4,2) 


- nbi5 


= 



(22) 

Again we have been careful to label the momenta where necessary. In fact, the colour factors 
of a given triangle pick up a minus sign (due to antisymmetry of the structure constants) under 
interchange of any of the three legs (2,3,5), as labelled in figure [3j BCJ duality (from eq. ([s])) then 
dictates that the corresponding numerators should obey similar relations. This implies that the 



first two terms of each relation in eqs. (21) and (22) cancel each other, leaving 



nbi = V i. (23) 
That is, the numerators for all bubble graphs should be zero. The numerators for the triangles are 



then fixed by the relations of eq. (20), each involving two boxes and a triangle. 



Having obtained the BCJ relations, one can show that these are automatically satisfied in the soft 
limit, in a Feynman gauge calculation. This is the subject of the following subsection. 



3.2 Soft limit 

In this section, we consider all possible infrared singular diagrams at one loop order, and show 
that their numerators explicitly satisfy the BCJ relations obtained in the previous section. IR 
singularities at one loop are generated by dressing tree level diagrams with soft gluon emissions. 
Considering once again the concrete example of 4-point scattering, we may write the full tree level 
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amplitude 



(24) 



where {n^;} is a suitable set of numerators (with colour factors {q} ) obeying the tree-level BCJ 
relations 

Us + Uu - Tit = 0, Cs + Cu - Ct = 0. (25) 

Note that we have absorbed coupling factors into the tree level numerators, so that all powers of the 



coupling in subsequent equations correspond to higher order corrections. Also, implicit in eq. (24) 



is the fact that the four-gluon vertex graph has been rewritten in terms of cubic graphs. The effect 
of dressing a tree-level amplitude Ax with a soft gluon exchanged between legs i and j is given by 



;Tj • Tj Xjj Axi 



where the eikonal integral factor 



Pi ■ Pj 



(26) 



as results from connecting two eikonal Feynman rules (eq. (11)) with a gluon propagator (note we 
use the Feynman gauge). The full infrared singular part of the one- loop amplitude can then be 
written |£] 



(1) 



IR 



T ■ T • 



A7 



(27) 



X i<j 

Here the first sum is over the three tree- level cubic topologies, and the second sum is over all pairs of 
external legs, where each pair is counted only once. There are no contributions from soft emissions 
which begin and end on the same external line, due to the fact that all outgoing hard particles 
are massless {pf = 0). Furthermore, there are no contributions from internal self-energy diagrams: 
such graphs contain additional hard propagators that remove the infrared singularity. 



We can now interpret eq. (27) from the viewpoint of BCJ duality. Firstly, the lack of internal and 



external self-energy graphs means that the numerators for these graphs may be set to zero in the 
soft limit. We may write this as 

hM = V i, (28) 

where the hat above the numerator implies that we are evaluating this only in the soft limit. Such 
numerators are defined only up to arbitrary terms containing at least one power of a soft gluon 
momentum k, as this kills off the infrared singularity generated by the integration over k. Equa- 



tion (28) indeed satisfies the requirement of eq. (23). 



^Note that we here adopt the phase conventions of eq. ([T]). Associating a factor i and ~i with each propagator 
and vertex respectively, one may absorb a further factor of 1 /i in the numerators to obtain an overall power of 
at L-loop order. Performing the double copy involves an extra numerator, and thus an additional explicit factor of i 
in the gravity amplitude. 

^Strictly speaking the eikonal integral of eq. (261 is zero, being a scaleless integral in dimensional regularisation. 
However, this is due to the cancellation of the IR divergence with a spurious UV pole which results from the 
replacement of quadratic propagators by linear ones, such that one may recover the IR divergence using an additional 
regularisation procedure. 



12 



Figure 5: Diagrams involving a soft gluon exchange between legs 3 and 4, where solid (curly) lines 
denote hard (soft) gluons respectively. 



Secondly, the terms appearing in eq. (27) can be grouped into sets of three, whose numerators 



satisfy the remaining BCJ relations. Interchanging the orders of summation in eq. (27), one may 
pick out a particular pair (i, j), to give 



Pi ■ Pj 



{27t)^ k'^ Pi • kpj ' k 



(29) 



where we have substituted in the explicit forms for Xj,- and Ax from eqs. (24) and (26). This is 



three separate terms, each of which can be interpreted as a soft limit of one of the graphs appearing 
in figure [3} T aking the pair (3,4) as an example, we may interpret eq. (29) as shown in figure [sj 
From eq. (|29[), we may associate with each term a collected numerator and colour factor 



One may then form the relations 



n, 



+ flu. 



nt. 



0, 



+ Cu 



,1-3 



Ct 



0, 



(30) 



(31) 



which follow from the fact that the tree level colour factors and numerators satisfy eq. (25), and 



that these are multiplied by a common factor in eq. (30). Equation (31) contains six different 



relations (the number of ways of choosing two external legs out of four), and correspond precisely 
to soft limits of the BCJ relations in eq. ( |20| ). Take, for example, the case shown in figure ([s]). 
These diagrams consist of soft limits of graphs (-B1), {B2) and (r3) in figure |3| and the relevant 
BCJ relation from eq. (31) thus corresponds to the second relation in eq. (20). 



We have now seen that all of the BCJ relations derived in the previous section (for the case of 
full QCD) are satisfied in the soft limit. For strict BCJ duality to hold, however, we must also 
demonstrate that the numerators are consistent with eq. ([5]), namely that appropriate sign changes 
occur under interchange of internal legs. That this is indeed the case is described in appendix [A} 



Armed with the above knowledge, we are now permitted to take the double copy of eq. (27) as 
prescribed in to give[^ 



M 



(1) _ _ -sp / f d^k (Pi-Pj)'^ rixUx 
IR - Z^Z^y^J J {27r)D k^ Pi- kpj-k X 



X Kj 



'^Here we have again adopted the phase conventions of eq. |ll - see the footnote on p. 12. 
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EE 



2 f d^k 2{pi-pjf 



i<j 



(27r)^ k"^ Pi ■ k pj ■ k x 



(32) 



If the double copy procedure is valid, this result should give the infrared singular parts of the 
4-point, 1 loop amplitude in GR. That this is indeed the case can be seen as follows. Firstly, the 
factor 



is, by the tree-level double copy procedure, a gravitational tree level graph for the s, t or u topol- 
ogy]^ This is then dressed by an eikonal integral, which corresponds to the appropriate gravitational 
generalisation of eq. (26). To check this, one may contract eikonal Feynman rules on legs i and j 



with the propagator of eq. (|18|) to get 

2 



T-grav. 



d^k p>'p'( 



PjPj 



{2tt)^ Pi ■ k \ Pj ■ k 

d^k 2{pi-pjf 
{2it)^ k"^ Pi ■ kpj ■ k^ 




D 



(33) 



where the minus sign in the second eikonal Feynman rule results from reversing the sign of the soft 
gluon momentum. Dressing the tree-level gravitational interaction with equation (33) is in exact 
agreement with eq. (32). 

Note that, as first observed in |4^, soft collinear singularities cancel after summing over all diagrams. 
Taking those which are associated with line i as an example, these are generated by the sum of 
eikonal integrals 

4.{pi-pjf "1 



E 



d^k 
- — -pr lim 
{2ti)D k^P, 



k'^ Pi ■ k Pj ■ k 



Applying momentum conservation 



Pi 



d^k 
\D 



lim 

(2^)^ fc^P. 



k^ Pi • k 



(34) 



(35) 



then gives zero on the right-hand side of eq. (|34|), up to non-singular terms. Note that the QCD 

d^k 



equivalent of eq. (34) (also including colour factors) is 

d^k 



Et T 

Colour conservation 



lim 



(27r)^ k^p, 



jPi-Pj) 
k"^ Pi ■ k Pj ■ k 



lim 



{2tt)D k^P, 



1 

k'^ Pi ■ k 



then gives 



Et. 



/ 



d^k 



lim 



(27r)'0 k^p, 



(Pi-Pj) 
k'^ Pi ■ k Pj ■ k 



Ci 



d^k 



lim 



(27r)^ k^p, 



1 



k"^ Pi ■ k 



(36) 
(37) 

(38) 



Recall that we absorbed gauge theory coupling factors into the tree-level numerators. Implicit in eq. (32 1 is that 
these have been replaced appropriately in the double copy to gravity. 
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where Ci is the quadratic Casimir associated with external hne i {Ci = Ca = Nc for pure glu- 
odynamics). This remains singular, and a simple physical way to interpret this is that collinear 
singularities can only depend on the quantum numbers of a single particle, which include the 
squared charge of the appropriate external line. In QCD, this is the quadratic Casimir operator in 
the relevant representation, whereas in gravity this is the squared 4-momentum, which vanishes if 
collinear singularities are to be present (i.e. if the leg is massless). The double copy procedure has 
here provided an interesting mechanism for the cancellation of soft-collinear singularities on the 
gravity side: the colour dependence in QCD gets replaced by additional momentum factors, which 
generate the necessary squared 4-momentum. Note that this means that there are singularities on 
the gauge theory side that vanish upon performing the double copy. We will see this happening 
more generally at two-loop order and beyond. 

We have now verified that the soft limit of GR is precisely reproduced upon taking the double 
copy of the soft limit of QCD at one-loop order. Although we here focused on the particular case 
of 4-point scattering, the argument is easily generalised to any number of external legs. Having 
examined the one-loop case, we proceed to two loops in the following section. 



4 Two loop analysis 

In the previous section, we showed that the known infrared singularities of GR are reproduced by 
double copying those of QCD at one-loop order. Application of the double copy relied upon the 
fact that the BCJ relations could be satisfied in the soft limit. In fact, solution of these relations 
at this order was rather straightforward, and relied ultimately on the fact that they are satisfied at 



tree level. Put another way, the soft BCJ relations of eq. (31 ) all had the form of a common eikonal 
factor (pi ■ pj) multiplying a tree-level relation. This simple structure will no longer be the case 
at two loop order. We will see that numerators for individual graphs, computed in the Feynman 
gauge, do not automatically satisfy the BCJ relations in the soft limit. However, this will turn out 
to be irrelevant for reproducing the known infrared singularities of GR. 



Our strategy will be the same as at one-loop order, and begins by writing down the BCJ relations 
in full QCD. Then one draws all possible soft topologies, and demonstrates that the resulting nu- 
merators satisfy the BCJ relations in the soft limit. Our experience at one loop tells us that we do 
not have to worry about most of the bubble diagrams. The numerators for these vanish as before, 
which will again be satisfied in the soft limit by the fact that internal and external self-energies 
are IR-suppressed. The only exception to this is the presence of self-energies associated with soft 
gluons, as we will see in what follows. There is a large number of diagrams at two-loop order and 
we do not collect them all here. Rather, we consider directly relevant soft topologies, and show 
that these can indeed be matched up with BCJ relations. 

We consider tree-level m-point scattering dressed by soft gluons up to two-loop order. Possible soft 
topologies can then link two, three or four external lines. Examples are shown in figure [6j where 
we label each topology. We may also write the eikonal integrals that go with each soft topology, as 

T.d.W = _ [ '^^^ / Pi ■ Pj Pk ■ Pi 

'^^^ J i2n)D J {2T:)DkH^pi-kpj-kpk-lpi-V 
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(Ik) 



(II3) 



(Y3) 







(II) (Y) (X) (O) 

Figure 6: Two-loop soft topologies in m-point scattering, where only four external lines are shown. 



where 



"ijk 
"ijk 

t(II) - - 

-^ij ~ 



f d^k r d^l 
J J2^J (2^ 

d^k 



Pi ■ Pj Pj ■ Pk 



r d^k f 
J WWJ 



r(Y) 

-^ij 



AX) 

^ij 



r(0) 



jPk_ 
(27r)^ 



_ f d^k f d^l 



k'^ Ppi ■ kpj ■ kpj ■ {k + I) Pk ■ V 

d^i Vt..p{Ki)ptp-^pi 

(27r)^ fc2 12 + ly p. . j^p. . ip^ .(k + iy 

_d^ {Pi-Pj? 

(27r)^ k^Ppi • kpi • {k + l)pj ■ kpk -(k + l)'- 
d^l V>..p(.k,l)p^p'^^p'; 
(27r)^ A;2 P {k + 0^ Pi ■ kpj ■ Ipj ■ k' 
{Pi-Pjf 



d^^l 



k'^ Ppi ■ kpi ■ {k + l)pj ■ Ipk ■ {k + iy 



d^k 

(2^7 {2Tr)i^ {kYlHl-k)^ Pi- kpj-k' 



V^,p{k,l)^0{k,l) 



(39) 



(40) 



is the three-gluon vertex, which in this case couples together only soft momenta. Our notation for 

the eikonal integrals specifies which indices are coupled together by soft gluon emissions, and note 

(Y) (Y) 

that the ordering of these indices can be important (e.g. I^j is not the same as Xj^ : the latter 
corresponds to a Y graph which is reflected with respect to the former) . The infrared singular part 
of the two-loop gauge theory amplitude can now be written as 



E 



xE{s,t,u} 



_{ijkl) 



^l-^ijkl ^ 



E(t. 

(ijk) 



^ k Hjk 



+ 



labc rpa rrib rrtc -T-C^s) \ 
/ ^i ^j ^k-^ijk ) 
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+ ^ ( (T, • Tjf 2^.1) + t5 T'j lip + T^^ T| iff^ + /"''^ / 



r(X) 



eabc jibdc rpa rpd -y-' 



nx 

X 



(41) 



where the notation {ij . . .k) denotes that one must sum over all distinct multiples, and we have used 
the colour vertex factor of eq. ([2]). This expression is obtained by dressing the tree- level interaction 
with all possible eikonal integrals, and is the two- loop generalisation of eq. (27). Furthermore, 
we have taken a factor of 
convention of eq. ([l]) . 



-1 out of the eikonal integrals, so as to make manifest the phase 



As in the one-loop case, we can match up terms in eq. (41) into sets of three, such that they 



correspond to the soft limit of a BCJ relation. There are two distinct classes of relation. Firstly, 
there are relations in which the same eikonal integral dresses different tree-level interaction graphs. 
This is the only scenario that was possible at one-loop order, and the relevant three graphs are 



obtained by picking a particular term in eq. (41 ), and keeping the sum over tree-level topologies x. 
As an example, consider the contribution 



(II) 

ij ' 



(42) 



obtained by selecting a particular pair in the first term in the second line of eq. (41). This corre- 



sponds to the graphs shown in figure [7][^ a). These indeed correspond to the soft limit of three full 
QCD graphs which enter a BCJ relation, namely those shown in figure [Tj^b). From eq. (42), one 
may associate a collected numerator and colour factor for each graph, according to 



.(II) 



n 



(II) 



(Pi -Pj) 



(43) 



where the hat once again reminds us that such numerators are defined in the soft limit, and may 
have arbitrary terms ~ 0{k,l) added. The numerators thus defined satisfy the relations 



nllD+nlll) 



ri- 



(II) 



0, 



(44) 



by virtue of the fact that this is satisfied by the tree-level numerators, and that these are all mul- 
tiplied by a common factor. 



It is clear that the above argument generalises to any set of three graphs which have the same soft 
topology Z. That is, one may form soft numerators 



n. 



j-iZ) 
-ij...k 



fix, 



where X) 



{Z) 



ij...k 



denotes the numerator of the relevant eikonal integral, such that 



nl^) + n. 



iz) _ AZ) 



0. 



(45) 



(46) 



This then corresponds to the soft limit of a BCJ relation. 
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(b) 





Figure 7: (a) Soft graphs corresponding to eq. (42); (b) full QCD graphs in the corresponding BCJ 
relation. 



The second class of soft BCJ relations at two loop order is more complicated, and consists of sets 
of three graphs in which the hard part of the interaction (in this case a single tree-level topology) 
is the same in each term, but the soft graphs are different. All possible examples are shown in 
figure [sj and the full (non-soft) BCJ triples can be obtained by replacing the soft gluon lines with 
solid lines, as in figure [Tj^a) and (b). Note that in figure [8 'c), the Y graph is related to itself, 
analagously to the case of the triangles considered in section 3.1[ For BCJ duality to be satisfied, 
there must exist soft numerators satisfying 



,(ll) 



n 



la) 
'x,ijk 



n. 



(O) 
x,ij 



corresponding to each of the sets of graphs in figure 



0, 
0, 
0, 

and where re 



(47) 



is the soft numerator 



associated with a diagram where soft graph Z dresses tree-level topology x. As is clear from 



the eikonal integrals of eq. (39), the numerators of the graphs in the Feynman gauge do not 



automatically satisfy these relations. The numerators are given by 



n 



3) 
x,ijk 

'^x,ijk 

(y) 

x,ij 

(II) 
x,ij 

~W 
'^x,ij 

■AO) 



re: 



re 



Pi -PjPj -Pknx, 
(u i\ 

iP 

{Pi ■Pjfria:, 



/iiyp{k,l)PiPjPjn^ 
■Pjfn^, 



^x ij = ^tiyp{^^ yyapik., I) K "-2 



(48) 



where we have used a tilde to denote the fact that these soft numerators do not respect BCJ duality. 
In order to find those that do, one must effect a generalised gauge transformation re — )• re. The 
general form of such a transformation is given in eq. ^ , and consists of redefining each numerator 
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Figure 8: Sets of three graphs, each corresponding to the soft hmit of a BCJ relation, such that 
the hard interaction is the same in each term. 
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through superpositions of denominator factors. If, however, all we care about is verifying the 
infrared singularities of gravity via the double copy, we do not have to solve the BCJ relations 
explicitly. To see this, note the numerator factors depend on soft momenta as 



n. 



(49) 



where K = (k, I) represents a soft momentum scale. Furthermore, all relevant denominator factors 
scale as at least 0{K). Thus, from eq. ([6| one must have 



n, 
fi 



3) 

x,ijk 

(II) 
x,zj 



n 



Is) 

x,ijk 



+ 0{K), 



n 



n. 



(II) 

X,lj 

■AX) 



X,l] 



and 



x,ijk ' 



''x,i3 ' 



x,tj 



0{K). 



(50) 



(51) 



We may summarise this more generally as follows. The BCJ dual numerators for soft graphs 
involving multiple single gluon emissions between pairs of external lines are the same as those in 
the Feynman gauge, up to corrections subleading in soft momentum. These corrections will not 
lead to additional IR singularities, so can be ignored in the soft limit. This is consistent with the 



to 



first two BCJ relations in eq. (47), which due to the subleading nature of n^^fji^ and iiX/^ amount 



n 



I3) 

x,ijk 



n 



II3) 

x,kji' 



n 



X,IJ 



n 



x,ij ' 



(52) 



as is indeed observed already in eq. ( 48 ) 



The BCJ dual numerators for graphs involving three gluon vertices off the eikonal lines are changed 
with respect to the Feynman gauge numerators, and are first order in soft momentum. Upon taking 
the double copy to gravity, these numerators will be squared, whereas the denominators are un- 
changed. Hence, by power counting, such graphs will not contribute infrared singularities in gravity. 
In other words, the very BCJ relations that one has to invest effort in solving are irrelevant for the 
double copy! This also tells us that there are infrared singularities in QCD that vanish when one 
takes the double copy to gravity. This "information loss" will, as we will see, be a feature at higher 
loop orders. 



The above remarks imply that we can take the double copy of eq. (41) by keeping only the graphs 



with multiple dipole emissions. Their numerators will be given by eq. (50), where we can safely ne- 



glect the corrections which are subleading in soft momenta. Performing the double copy procedure 
gives a 2-loop gravity amplitude |^ 



M 



(2) 
IR 



I) E 

x&{s,t,u} 



E7-(l|4),grav. v;^ .^(||3),grav. v;^ /^(||),grav. 
-^ijkl -^ijk 



I 7-(^).gi'av. 



(ijkl) 



{ijk) 



(ij) 



Tlx f^x 



X 



(53) 



For the overall phase factor, see the footnote on p. 12. 
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where 



•7-(l|4),grav. 



ijk 



T-dD.grav. 

j-{X),grav. 
"ij 



d^k 
(2^ 

d^k 
(2^ 



d^l ^{Pi ■ Pjf {Pk ■ Pi? 
(27r)^ k"^ P Pi ■ k pj ■ k pj. ■ I pi ■ V 

d^l ^{Pi-Pj?{Pj -Pk? 

(27r)^ A:2 ^2 p. . p^. . kpy {k + I) pj, ■ V 

d^l 4(K-pj)4 

(27r)^ A:2/2 

Pi ■ kpi ■ {k + l)pj ■ kpj • {k + I) 

_d^ 4.{m-Pj)^ 

(27r)^ k^ Ppi-kpi-{k + I) pj -Ipj-ik + l)' 



(54) 



This indeed agrees with an exphcit calculation using the known gravitational eikonal factors, ob- 
tained using the Feynman rules of eqs. (17) and (18). We have thus shown that, at two loop order, 
the infrared singularities of GR are consistent with those of QCD via the double copy procedure. 



We can write eq. (54) in a more recognisable form as follows. Firstly, one has 



j(l|4),grav. 



7-grav.^grav. /c-c-x 
-ijkl —-^ij -^kl ' y^^) 

where the right-hand side contains the product of two gravitational one-loop eikonal factors from 
eq. ( 33 ) . By collecting terms in eq. ( 54 ) , we can write the entire right-hand side in terms of one- loop 

^Pi ■ Pj)^ {Pj ■ Pkf 



integrals. In particular, one has 



j(-f-f3),grav. j-(Il3),gTav. 
ijk ' kji 



d^k 



d^l 



{2tt)D J {27r)D k^Ppi-kpk-l Pj ■{k + l) 

d^'k f d^l ^{p^■Pjf{pj■Pk? 
{2ti)D J {2n)D k^Ppi • kpk ■ Ipj ■ kpj -V 



1 1 



where in the second line we have used the eikonal identity 



1 



Pj -{k + l) 



1 1 

+ 



Pj ■ k Pj ■ I 



1 



One thus has 



Also, one has 



j(ll):grav. 2-('''^)'gi'^*^- 
ij "T ij 



j(^^3),grav. j{IIz),gTmr. 
ijk ' kji 



Pj ■ kpj ■ I 
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ij jk • 
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jPl ^{Pi-Pj)^ 
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jPl KPi-PjY 

(27r)^ A;2 P p^ ■ kpi ■ {k + l)pj ■ kpj ■ I 
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d^l ^{Pi-Pjf 



1 1 



1 1 

p^^k^pT-n 



(56) 

(57) 
(58) 
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2\ ij , 



(59) 



where we have again used the eikonal identity of eq. (57), as well as relabelling k ^ I. Putting 



things together, we may rewrite eq. (53) using the results of eqs. (55, 58, 59) to give 



M 



(2) 
IR 



2 



x(i{s,t,u\ 



En 

(ijkl) 



grav..T-grav. 
ij -^kl 



■grav. 
ij 



(60) 



{ijk) (ij) 

where the notation (ijk) denotes summing over all triples, such that the ordering of i and k is 



unimportant. The total soft prefactor (contents of the square bracket in eq. (60) is easily checked 
to be the second order term in the expansion of 



exp 



(ij) 



(61) 



as expected from the known exponentiation and one-loop exactness of gravitational infrared diver- 
gences 49,50 52 . This completes our analysis of BCJ duality and the double copy in the soft limit 



at two-loop order. In the next section, we generalise our remarks to all loop orders. 



5 General remarks 

In the previous two sections, we have seen that the structure of infrared singularities in QCD 
matches on to those of GR after applying the double copy procedure. At one loop, we were able to 
apply the double copy due to the fact that the Feynman gauge numerators for soft graphs automat- 
ically satisfied the appropriate BCJ relations in the soft limit. At two loops, BCJ relations could be 
separated into two classes: (a) those involving graphs whose soft topology was the same, but whose 
underlying hard tree- level topology was different; (b) those involving graphs sharing the same hard 
interaction, but having different soft topologies. BCJ relations of class (a) were still automatically 
satisfied, due to the fact that the underlying tree-level numerators satisfy BCJ duality. Relations of 
class (b) were more complicated. Numerators for graphs with no three gluon vertices off the exter- 
nal lines (which we refer to as dipole-like graphs in the following) could be taken to be the same as 
the Feynman gauge results up to subleading corrections in soft momenta, which are irrelevant from 
the point of view of infrared singularities. Numerators of graphs involving three- gluon vertices off 
the external lines did not automatically satisfy BCJ relations, and thus would have to be modified 
by generalised gauge transformations in order to write down a BCJ-dual representation of a QCD 
amplitude in the soft limit. However, these graphs are at least linear in soft momenta, and thus 
vanish upon taking the double copy to gravity. This allowed us to verify that the IR singularities 
of gravity are correctly reproduced by double copying the QCD results at two-loop order, without 
having to explicitly solve the BCJ relations. 



The aim of this section is to argue that this argument generalises to all loop orders. This is possible 
because we have already seen all of the necessary ingredients at two loop order. Firstly, the fact 
that the infrared limit BCJ relations fall into the two classes given above is generally true, indepen- 



dent of the loop order " Then BCJ relations of class (a), and involving a given soft topology (Z), 



At one loop, as we have seen, only class (a) occurs. 
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take the generic form of eq. (46), with numerators given by eq. (45). As at two loop order, these 
relations are satisfied by virtue of the fact that the tree-level numerators satisfy the BCJ relation, 
and are multiplied by a common factor. 



BCJ relations of class (b) are again more complicated, and can be split into two further subclasses: 
(i) those involving two dipole-like graphs and a graph containing at least one three gluon vertex 
off the external lines; (ii) those involving three graphs with at least one three gluon vertex off the 
external lines. By the same power counting arguments that were used in the previous section, the 
numerators of dipole-like graphs are the same as their Feynman gauge counterparts up to correc- 
tions which are subleading in soft momenta (i.e. which do not contribute infrared singularities). 
Two dipole graphs which enter the same BCJ relation are related by a permutation of two gluon 
emissions on an external line (e.g. figure [8](a)) Thus, BCJ relations of subclass (i) set the 
numerators of such graphs to be equal, which is indeed satisfied in the Feynman gauge, which asso- 
ciates with a single soft gluon emission between lines i and j a contribution 2{pi -pj), independently 
of any other gluon emissions. BCJ relations of subclass (ii) are not satisfied by the Feynman gauge 
numerators for the relevant graphs. However, these numerators are all at least 0{K) (where K is 
an arbitrary soft momentum), and remain so after performing a generalised gauge transformation 
in line with eq. ([6|. Thus, these graphs do not give infrared singularities after performing the 
double copy. One thus does not need to solve the BCJ relations explicitly in order to generate the 
IR divergences of the gravity amplitude. From three-loop order in the Feynman gauge, one must 
also consider graphs involving the four-gluon vertex off the eikonal lines. This also gives rise to 
numerators which involve non-zero powers of soft momenta, after rewriting such graphs in terms 
of cubic topologies, according to the usual BCJ procedure. 



The above remarks allow us to generalise eq. (60) to an arbitrary loop order, as 
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•jgrav. 
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(62) 



In this formula. If •' is the eikonal integral factor due to the sum of all dipole emissions that 



link lines 11,12 ■■ .im- The second sum in the square brackets in eq. (62 ) is then over all multipoles 



{ii ■ ■ ■ i-m) that are linked by such dipole emissions. The first sum is then over all possible numbers 
m of external lines. As in the two-loop analysis of the previous section, one may collect terms in 



eq. (62) into products of one- loop eikonal integrals, via multiple applications of the (higher-order) 

2pi ■ Pj 



eikonal identity, and thus rewrite eq. ( 62 ) as 491 
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(63) 



The tree-level amplitude dressed by soft gravitons to all orders, and obtained via the double copy, 
is then given by 
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IR 
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IR 



exp 
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iij) 
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d k 2pi • Pj 
(27r)^ k^Pi -kpj-k 



(64) 



In the language of |29l 58 59 , such diagrams are in the same multiaprton web. 



^^Here we have absorbed factors of i into the eikonal integrals, rather than show these explicitly as in eq. 
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in agreement with the known all-order structure of IR divergences in GR 49 . Note that, once we 



had reached eq. (62), the final result was guaranteed: what mattered was that the double copy 
procedure correctly reproduces the fact that only dipole-like graphs appear in gravity. The sum 
over all such graphs automatically leads to the exponentiation of the one-loop corrections 
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This completes our argument that the all-order structure of IR divergences in gravity is consistent 
with IR singularities in QCD, via the application of the double copy procedure. This can therefore 
be taken as all-order evidence for the double copy conjecture, albeit in a particular limit. Some 
further comments are in order. Firstly, the general argument at 0{a^) exhibits a property already 
remarked upon at two loops, namely that many singularities cancel upon taking the double copy 
to gravity. This means that the gravity side of the correspondence cannot be used to constrain 
singularities on the gauge theory side. For this reason, the question posed earlier in the paper 
regarding whether the QCD dipole formula has a gravitational origin appears to have a negative 
answer: the very singularities which could occur as corrections to the dipole formula vanish when 
we move to gravity. 

Secondly, we have here focused on the case of pure Yang-Mills theory and General Relativity. How- 
ever, the argument of this paper can also in principle be applied in supersymmetric gauge theories 
/ supergravity. Here we remark that reference [l6| obtained amplitudes in > 4 supergravity 
using the double copy procedure applied to gauge theory amplitudes in = 4 Super- Yang- Mills 
theory coupled with < < 4 Super- Yang-Mills theory. One check on this calculation was the 
demonstration that infrared singularities, after the double copy, were consistent with exponentia- 
tion on the gravity side up to two loop order. The results of this paper generalise this to all loop 
orders, and non-supersymmetric theories. In supersymmetric theories, it may be possible to extend 
our arguments beyond the pure soft limit, as non-trivial information can in principle be obtained 
from infrared singularities |66[|67| . 



6 Conclusion 



In this paper, we have examined the soft limits of QCD (strictly speaking, pure gluodynamics) 
and GR, and showed that infrared singular contributions to amplitudes in both theories match up 
with each other upon using the double copy procedure of [2||3] . The structure of IR divergences in 
QCD scattering amplitudes is still subject to some uncertainty, and our current state of knowledge 

By 



can be expressed in terms of the dipole formula of 31-34 , plus possible corrections 135 -39 



contrast, the structure of IR singularities is known exactly in GR 49-52 
in terms of one-loop graphs only. 



where they exponentiate 



Being able to take the double copy relies on the gauge theory amplitudes displaying manifest BCJ 
duality [s]. At one- loop order, we saw that amplitudes in the soft limit, as calculated in the Feyn- 
man gauge, automatically satisfied BCJ duality. This was not the case at two-loop order. By 
power counting, however, this turns out to be irrelevant for graphs consisting of multiple dipole 
emissions. These are the only graphs that survive upon taking the double copy, so that one does 



^^This is the same argument that occurs in the exponentiation of one-loop soft corrections in QED 
theory, however, one-loop exactness is broken due to the presence of fermion bubbles. 



^0||49 



In that 
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not need to solve the BCJ relations explicitly. We could thus use the double copy procedure to 
"predict" the structure of IR singularities in GR, finding exact agreement with the known results. 
Our arguments imply that many singularities vanish upon taking the double copy. If this were not 
the case, one could have used singularities in gravity to constrain possible corrections to the QCD 
dipole formula, or even to provide an underlying gravitational explanation for this. Nevertheless, 
the arguments presented in this paper constitute all-loop level evidence for the validity of the double 
copy conjecture, which may be more significant in supersymmetric contexts. 

Note that we have only considered the case where the hard interaction consists of tree-level scat- 
tering. In principle, this could contain higher loop orders. To show that BCJ duality is satisfied 
by the full amplitude then requires BCJ-dual numerators for the hard interaction, which is not 
possible to all orders without a full proof of BCJ duality in QCD. 



There are a number of further questions that can be addressed. It may be possible, for example, 
to reinterpret our results using the manifestly BCJ-dual effective Lagrangian of It would also 
be interesting to examine the full implications of the present analysis in supersymmetric contexts. 
Finally, a thorough investigation of the role of BCJ duality in non-supersymmetric gauge theory 
away from the soft limit would potentially provide new insights into QCD and / or gravity. An 
intermediate step in this regard might be to extend the present analysis to next-to-eikonal order, 
using the technology of (51 
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A Numerator reflection properties under particle interchange 

In this appendix, we show that the numerators obtained for the 1-loop triangle graphs in the soft 
limit obey the BCJ reflection property of eq. ([s]), namely that if the colour factor changes by a sign 
under interchanging two (external or internal) legs, the numerator also does so. It will be sufficient 
to consider a particular example, for which we choose the graph labelled (Tl) in figure [Sj 

In the soft limit, this graph may be redrawn as shown in figure [9j where we have assigned a 
particular momentum flow. The labels on the legs can be translated into physical momenta as 
follows: 

3 = k, l = 2=pi, A = 5 = p2, (65) 
using the fact that the soft gluon momentum /c — )• 0. The numerator for this graph in the soft limit 



is given by eq. ( 30 ) as 

nTi(l, 2, 3, 4, 5) = (pi • P2)ns{2, 5), (66) 

where in the tree-level numerator for the s-channel graph we note explicitly that this has the mo- 
menta 2 and 5 as its incoming legs. 
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Figure 9: The soft limft of the triangle diagram (Tl) in figure [sj 



The colour factor of this graph picks up a minus sign under the interchanges 2 -f-)- 5, 2 -f-)- 3 and 
3 -f-)- 5. Making the first interchange and imposing momentum conservation gives 

3 = k, 1 = 5= P2, 2 = A = pi. (67) 

The numerator factor itself maps to 

nTi(l,2,3,4,5) ^ (pi-p2)n,(5,2). (68) 

We may now use the fact that the tree-level numerator is BCJ dual, and thus must satisfy 

n,(2,5) = -n,(5,2). (69) 



Substituting this into eq. ( 68 ) then yields 



nri(l, 2, 3, 4, 5) ^ -nTi(l, 2, 3, 4, 5) 



as required. 



Next, we consider the interchange of legs 2 3, which modifies the graph shown in figure |9] to give 
that of figure [Toj^a), where we have included relevant physical momenta, after imposing momentum 
conservation. The resulting graph is not infrared singular as — )• 0. Indeed, the k dependence of 
the original graph 

d'^k 1 



{2ttY {Pi - kY {p2 + kf 



have been shifted to 

d^k 1 



(2'KYk'^ {k+piY{k-p2f' 

The infrared singular region is found by shifting k — )• k + pi, after which one obtains the graph 
shown in figure 10 b), which has k dependence 

d-^k 1 



(2^)'^ A;2 {k + piY {k - P2Y' 



This corresponds to the original graph, with pi ^ P2- Put another way, in order to make figure 10 ^b) 
have the same momentum fiow as figure [9j we must fiip the two incoming lines. Thus, one has 

nTi(l,3,2,4,5) =pi •p2n,(5,2), (70) 
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(a) (b) 

Figure 10: (a) Diagram obtained by interchanging legs 2 o 3 in figure [oj (b) Diagram obtained by 
shifting — )• k + pi. 



where we have reinstated the notation used in figure [9j Again using eq. ( 69 ) , we find 

nTi(l, 3, 2, 4, 5) = -nTi(l, 2, 3, 4, 5) (71) 
as required. The case of 3 o 5 is directly analagous. 

Here we have considered a triangle built upon the s-channel tree-level graph. A similar analysis 
can be used for the t- and u-channel cases, as well as triangular subgraphs at higher loop orders. 
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